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Abstract 

The two-point correlation functions of statistical models show in 
general both poles and cuts in momentum space. The former corre- 
spond to the spectrum of massive excitations of the model, while the 
latter originate from interaction effects, namely creation and annihi- 
lation of virtual pairs of excitations. We discuss the effect of such 
interactions on the long distance behavior of correlation functions in 
configuration space, focusing on certain time-slice operators which are 
commonly used to extract the spectrum. For the 3D Ising model in 
the scaling region of the broken-symmetry phase, a one-loop calcula- 
tion shows that the interaction effects on time-slice correlations is non 
negligible for distances up to a few times the correlation length, and 
should therefore be taken into account when analysing Monte Carlo 
data. 
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1 Introduction 



When studying a statistical model, one is often interested in determining the 
spectrum of massive excitations, i.e. the eigenvalues of the transfer matrix. 
For many interesting models this cannot be done exactly, and one has to 
rely on approximate methods or numerical Monte Carlo calculations. 

The observables which are most suitable to investigate the spectrum of 
a model are the two-point correlation functions of operators: their long 
distance behavior is directly related to the spectrum. This is especially 
evident in momentum space, where each pole of the correlation function 
corresponds to a massive excitation, and therefore to an eigenvalue of the 
transfer matrix. 

In general, however, the momentum space correlators will have not only 
poles but also cuts, signaling the possibility of creating and annihilating 
virtual pairs of excitations. For example if the effective Hamiltonian for the 
order parameter <f> includes a c/> 3 interaction, the Feynman diagram 



will generate a cut in the {(j)(j>) correlator in momentum space. 

Since, in general, Monte Carlo simulations give direct access to config- 
uration space correlation functions rather than their Fourier transforms, it 
is interesting to study the effect of such interactions on the long distance 
behavior of configuration space correlators. The purpose of this work is 
to compute this effect in the scaling region of the broken symmetry phase 
of the 3D Ising model, where actual calculations can be performed in the 
framework of renormalized Euclidean quantum field theory. 

A strong motivation for this analysis is provided by the continuous im- 
provements in the accuracy of Monte Carlo simulations: recent advances in 
both computer performances and simulation algorithms allow us to obtain 
numerical data of unprecedented precision. Their analysis requires more 
sophisticated theoretical tools as finer effects become observable. We will 
find that the effects that are the object of this study typically account for 
about one percent of the correlators in the region of physical interest: this 
is actually an order of magnitude larger than the statistical uncertainties 
typical of recent Monte Carlo studies of the 3D Ising model. 

The paper is organized as follows: in Sect. 2 we introduce time-slice 
operators, which are particularly suitable for the Monte Carlo study of the 
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spectrum of a statistical model. In Sect. 3 we compute the correlators of 
such operators in the Ising model using (p 4 field theory at one loop in three 
Euclidean dimensions. In particular the interaction effects can be evaluated 
and expressed in terms of exponential integral functions. In Sect. 4 we 
comment on the relevance of these effects to certain universal amplitude 
ratios, while Sect. 5 is devoted to the discussion of the results and their 
implications for the analysis of Monte Carlo data. 



2 Time— slice correlators 

The spectrum of massive excitations of a statistical model can be obtained 
by studying the long distance behavior of two-point correlation functions. 
Consider for example the Gaussian model, with Hamiltonian 

H = f ^x^-d^d^+^ct^ . (1) 

The two-point correlation function is 

f d d p e^-y) 
{<j){x)<j)(y)) = / 77^-^^— — 2 

1 /\ X - yW 1 - 4 / 2 

where Ki_ d / 2 is a modified Bessel function. Therefore asymptotically for 
\x — y\ — > oo 

(4>(x)4>(y)) ~ const \ x - e -™\*-v\ . (3) 

We see that the long distance behavior of correlators can be used to extract 
the value of the mass m, which in this case is obviously the only state in the 
spectrum. 

In practice, it is more convenient to define time-slice operators 

S (tf = id~T J dxi...dx d -i(j>(xi,...,x d -i,t) (4) 

to obtain a purely exponential behavior of correlations. For example in the 
Gaussian model it is easy to see that 

(S(t)S(0)) = j^e-^ . (5) 
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When considering a non-trivial model, it is customary to generalize Eq.(||) 
to 

(S(i)S(O)} = $> fe e-™*l*l . (6) 

k 

By fitting the values of the time-slice correlations with Eq.(^) one can ex- 
tract the values of a certain number of low-lying states, depending on the 
precision of the available data. 

Each exponential in the r.h.s. of Eq.(||) corresponds to a pole in the 
Fourier transform of {(p(x)4>(y)) . However, since a non-trivial model cer- 
tainly involves interactions, Eq. must be modified to take into account 
their contribution: in the next Section we will compute the time-slice cor- 
relator (S(0)S(t)) for the 3-D Ising model in the broken symmetry phase, 
where a </> 3 interaction is present in the effective Hamiltonian so that cuts 
appear in the two-point function already at one loop. The calculation will 
teach us how to modify Eq.(^) to take into account the effect of produc- 
tion/annihilation of virtual pairs of excitations. 



3 The case of the 3D Ising model 

It is a widely accepted conjecture that the 3D Ising model is in the same 
universality class as (j) 4 field theory. This allows us to use renormalized 
3D quantum field theory to study the Ising model in the scaling region, 
where lattice effects become negligible and universality holds. This program 
was initiated by Parisi [|| and has been vigorously pursued to study several 
aspects of the Ising model ^, [|. The agreement between field theoretical 
calculations and Monte Carlo results is satisfactory. 

Therefore, from now on we will consider the 3-D Euclidean field theory 
defined by the action (effective hamiltonian, in the language of statistical 
mechanics) : 

We are interested in the two-point connected function 



d 6 x 



(7) 



G(x -y) = (<f>(x)<f>(y)) - (<f>(x)) (<j>(y)) . (8) 

The perturbative expansion must be performed around one of the stable 
classical solution, say <j> = v. Defining a fluctuation field ip = <f> — v a if 3 
term appears in the Lagrangian, with a coupling proportional to yj~g (for 
details about the perturbative expansion in the broken symmetry phase see 
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e.g. Ref. j3|). The correlation function (||) is then given at one loop by the 
sum of the following Feynman diagrams: 



G(x-y) 



+ 



n 



+ 



+ 



(a) (b) (c) 

Using dimensional regularization we find 



(d) 



Gd(x - y) 



d d p 
(2tt)° 



G d (p)e 



ip(x-y) 



(9) 



with 



G d (p) 



p 2 + m 2 



+ 




d/2 



r(i 



(p 2 + m 2 ) 2 



+ 



3m 2 g 



F d {p) 



(p 2 + m 2 ) 2 



and 



F d (p) 



r(2 - d/2) 
(4vr) rf / 2 



+ x(l — x)p z 



d/2 -2 



(10) 



The first term in the r.h.s. of Eq.(|l0|) is the tree level contribution 



with a pole in p 



-m 



The remaining diagrams are interaction effects: 



the second term, corresponding to diagrams (b)+(c), produces shift in the 
location of the pole, i.e. a quantum correction to the physical mass. The 
third term is the interesting one for our purpose: besides providing another 
quantum correction to the physical mass, it has a cut in p 2 = —Am 2 . Indeed, 
this term corresponds to diagram (d) in the expansion of G(x — y) namely 
to the production and annihilation of a virtual pair of particles. 

The theory must be renormalized to be compared with experimental or 
Monte Carlo results. Notice that analytic continuation to d = 3 gives a finite 
expression for G{p) without need for any subtraction. This is a peculiarity 
of dimensional regularization for odd d which disappears when higher loop 
effects are taken into account. It implies that the renormalized parameters 
in the MS scheme coincide at one loop with the bare parameters. 

However the MS scheme is not particularly suited for direct compari- 
son of field theoretic calculations with experimental or Monte Carlo data. A 
more convenient scheme, where the renormalized parameters of the field the- 
ory have a direct lattice interpretation was introduced in Ref. || . The renor- 
malized parameters <j)R, run, gn in this scheme were computed at three- loop 
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order in [Q], to which we refer the reader for the definitions and expressions 
of the renormalized parameters. In particular the renormalized mass m R is 
defined as the momentum space correlation at p = 0, and coincides with the 
inverse second-moment correlation length. Here we just need the expression 
of the connected two-point function of renormalized fields: 

G R (x -y) = (M*)My))c = I ^Gfl(p)e*^ (12) 
with, denoting with ur = Qr/tur the dimensionless renormalized coupling, 
p , \ ur \ I 1 m R u R 1 



3m R UR 1 1 



„2 



8. (^ + m|)W arCtan l^Ji (13) 
Defining the time-slice operators as 

S(t) = j2 J dxxdx2<j> R {xx,X2,t) (14) 

we have 

(S(t)S(0)) c = ±J ^e^G fl (0,0,p) (15) 
After a simple calculation we obtain for t > 



(*(t)*(0)> = 



1+ l^( 24l0g3 " 27 ) 



TTH. / Z2 ■ ( 16 ) 

l6nL z m R J2m R A _ 



■r 



where m p h is the physical mass, defined as the location of the zero of the 
inverse correlator in momentum space G _1 (p) 



2 2 

m ph = m R 



1 + S (13 - 12l0g3) 



(17) 



The integral appearing in Eq.(p^) can be expressed in terms of exponential 
integral functions to give [J] 

1 The appearance of an exponential term e~ 2mt in Eq.([l8]) could be misleading: due to 
a cancellation between this term and the Ei functions the asymptotic behavior of the sum 
in square brackets is actually e~ 2mt /t. 
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{S(t)S(0))c = 
167rL 2 m/j 



1 



2m R L 2 



-m ph t 



1 + 



1287T 



(24 log 3 -27) 



g-2m ph i 4 + 2 

___ + e -^*Ei(-Tnp fc t) 



2 ™^V ^Ei(-3m p ^) - Ei(-2m p/l i) 



(18) 



Eq.(18) is our main result: it gives the contribution of interaction effects 
to the correlation function of time-slice operators as a modification of the 
simple exponential behavior (0). 



4 The universal amplitude ratio £/£2nd 

Universal amplitude ratios are certain dimensionless combinations of observ- 
ables that are predicted to be universal at criticality (for a comprehensive 
review see Ref . || ) . Among these it is of particular relevance to us the ratio 
£,/£,2nd of the second-moment correlation length and the "true" correlation 
length (i.e. the inverse of the physical mass). This ratio defines two uni- 
versal amplitude ratios, F+(F-), when the critical limit is taken from the 
symmetric (broken symmetry) phase. It can be shown Q that the presence 
of higher masses in the spectrum implies F± > 1. However the correspond- 
ing analysis in 4 theory shows that the converse is not true: a value of F- 
greater than one does not necessarily indicate the presence of higher mass 
states, but can simply be a signal of non-trivial interaction effects like the 
ones studied in this work. 

In fact, since the renormalized mass is defined as lf&nd (see ||), 
from Eq.(|l7|) we see that at one loop in </> theory, || 

F_ = 1 - y||^ 13 - 12 log 3) = 1.00668(3) , (19) 

in good agreement with the Monte Carlo result H 

F_ = 1.009(5) (20) 

Moreover, it is easy to see that the non-trivial contribution to F^ comes 
exclusively from diagram (d) in the expression of G(x — y), namely the 
diagram which produces the cut. 
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It would be interesting to study similar effects in the symmetric phase: 
however we expect them to appear with the diagram with the diagram 




and therefore only at two-loop level. This provides a qualitative explanation 
for the fact that in the symmetric phase the corresponding amplitude ratio 
F + is known to be much smaller than F_ : a strong coupling expansion gives 
1 

F + = 1.00023(5) (21) 
while Monte Carlo calculations give an upper bound |?| 

F + < 1.0006 (22) 



5 Discussion 



The relevance of the effect we have just computed can be best appreciated 
by considering the ratios 

(5(t + l)5(0)) c 

^ ) = - log msWc (23) 

For a purely exponential behavior, R(t) is identically equal to m p h, while 
R(t) as given by Eq.(|i~8|) is plotted in Fig. 1, where we have used the Monte 
Carlo estimate 

u* R = 14.3(1) (24) 

for the value of the dimensionless renormalized coupling in the continuum 
limit and we have set m p h = 1 {%■&■ we are measuring distances in units 
of the correlation length). 

The figure shows that R(t) is appreciably different from 1 for distances 
of the order of a few times the correlation length. For example at t = 1 the 
interaction effect is ~ 0.8% of R. State-of-the-art Monte Carlo simulations 
give statistical uncertainties about ten times smaller for the same quantity 
jjj. The magnitude of the effect becomes comparable with the statistical 
uncertainties at t ~ 2.5. 



2 The variation of ur is anyway very slow in the whole scaling region, being governed 
by corrections to scaling. 
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Figure 1: The ratio R(t) defined in Eq.(2S) as predicted by Eg. fllty for 
m p h = 1. A purely exponential behavior would give R(t) = 1 identically. 



The specific form of Eq.fllq) suggests one more reason why the correction 
must be taken into account when analysing numerical data. In fact, the be- 
havior of the correlators when interaction contributions are included mimics 
very closely the contribution of a higher mass in the spectrum. Consider 



the ratio fl23|) when the contribution of two poles is included and the cut is 
neglected: 



R(t) = -log 



-m(t+l) 



+ ae' 



-m'(t+l) 



-rat 



+ ae' 



-ra't 



(25) 



We have verified that by adjusting the parameters a and m'/m one can make 
R(t) look very similar to the cut contribution: relative uncertainties in the 
data of less than one part in 10 4 would be needed to resolve the difference. 
The "best fit" value is about m'/m = 2.4. Therefore the effect of interaction 
can easily be mistaken for a higher mass state with m! ~ 2.4 m. 

A more complete investigation of these issues, including a high-precision 
Monte Carlo analysis, is currently being pursued and will be published else- 
where lie 
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